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System and Method for Shaping Ultra Wide Bandwidth Signal Spectrum 
Field of the Invention 

[01 ] This invention relates generally to wireless communications, and more 
particularly to shaping the spectrum of transmitted signals in ultra wide 
bandwidth communication systems. 

Cross Reference to Related Application 

[02] A claim of priority is made to U.S. provisional application serial no. 
60/451,577, "UWB Communication system with shaped signal spectrum by 
pulse combination," filed March 3, 2003. 

Background of the Invention 

[03] With the release of the "First Report and Order," Feb. 14th, 2002, by the 
Federal Communications Commission (FCC), interest in ultra wide bandwidth 
(UWB) communication systems has increased. The IEEE 802.15 standards 
organization, which is responsible for personal area networks (PAN), has 
established a task group, TG3a, to standardize a high-data-rate physical layer 
based on UWB. One of the most important requirements is compliance with the 
FCC spectral mask for indoor channels shown in Figure 1 . The spectral mask 
shows power limits as a function of frequency. 

[04] In addition to this requirement, it is also required by IEEE 802. 1 5 that 
UWB systems do not interfere with existing wireless systems, such as 802.15.1 
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(Bluetooth), 802.15.3 (personal area networks), 802.15.4 (Zigbee) and 802.1 la 
and 802.1 lb (wireless local area networks). 

[05] Furthermore, UWB systems should also be resilient against interference 
from these existing wireless systems. In addition, UWB systems should also 
resist interference from home appliances, such as microwave ovens, and other 
narrowband interferers. 

[06] All these requirements impose additional constraints on the spectrum 
shaping. Minimizing interference, in both receivers and transmitters, requires 
placing minima on the poser in the transmission spectrum. Rigorously speaking, 
the spectral mask demands very low power in certain frequency bands, as well 
as matching filters. A null in the transmit spectrum also implies a null in the 
transfer function of the matching filter, and thus interference suppression at this 
frequency. 

[07] While some interference, e.g., microwave radiation, is at fixed 
frequencies, other interference is at variable frequencies, e.g., the different 
bands of 802. 1 la, or at frequencies that cannot be predicted a priori. Thus, it is 
necessary to be able to shape the spectrum adaptively. 

[08] For impulse radio signals, pulse-position modulation, (PPM) and pulse 
amplitude modulation (PAM) are the most popular signaling techniques. These 
techniques are combined with time hopping (TH) for multiple-access format. In 
TH, each symbol is represented by a series of "basis pulses." The location or 
delay of the pulses is determined by a time-hopping code. The time-hopping 
code is generated as a pseudo-random sequence. 
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[09] In general, the spectrum with random code generation can be quite 
different from the FCC mask. As a resuh, in order to be compHant with FCC 
regulation, with a scaling of transmission power, the maximum allowed 
transmission power can be significantly restricted. Therefore, any signal 
shaping method that uses a basic, single short pulse must explicitly take the 
FCC mask into consideration. 

[010] Under certain conditions the spectrum of the transmit signal becomes 
identical to the spectrum of the basis pulse. Therefore, there is a need to modify 
the spectrum of the basis pulse to fit the requirements of both the FCC mask 
and industry standards. 

Summary of the Invention 

[Oil] The invention shapes the spectrum of a transmitted signal in an ultra wide 
bandwidtii (UWB) system so that the transmitted signal meets regulatory 
constraints, as well as practical requirements of interference suppression. 

[012] The system according to the invention uses a linear combination of 
weighted, time-shifted (delayed) basis pulses in order to achieve a desired 
spectrum. Characteristically, the system uses a bank of basic pulses, preferably 
with different time-frequency duration, to approximate a given mask of general 
shape with a small number of parameters. 
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[013] The weights and delays are selected according to criteria related to 
spectral efficiency, adherence to the FCC spectral mask, and other restrictions 
like maximum number of rake fingers. 

[014] The method for determining the weights is based on a non-linear search, 
initialized by the solution of a quadratic optimization problem. Two ways of 
approximating the original problem with a quadratic formulation are described. 
Additionally, the invention provides combinatorial optimization to perform 
discrete optimization over pulse positions (delays), as well as off-line selection 
from a basic pulse set. 

Brief Description of the Drawings 

[015] Figure 1 is a graph of a spectral mask used by the invention; 

[016] Figure 2 is a block diagram of a system for shaping an ultra wide 
bandwidth spectrum according to the invention; 

[017] Figure 3 is a block diagram of the system of Figure 1 further including a 
set of oscillators according to the invention; 

[018] Figure 4 is a block diagram of a feed-forward neural network used by the 
invention; and 

[019] Figure 5 is a block diagram of feed-backward neural network for 
adaptive optimization with back-propagation according to the invention. 
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Detailed Description of the Preferred Embodiment 

[020] As shown in Figure 2, our invention uses a linear combination 1 10 of a 
set of basis pulses 101 for shaping of a spectrum a transmitted impulse radio 
signal, s(t) 102, for example, a transmitted ultra wide bandwidth (UWB) signal. 
The basis pulses 101 are generated using circuits so that they are spread over 
the frequency spectmm with pseudo random delays, 

[021] Prior to the combining 1 10, a set of optimized filter and delay lines 120, 
hi(t), / = 1 , . . . , Z, is applied to the pulses to weight and delay the generated basis 
pulses to conform to a predetermine spectral mask. 

[022] The combination of the filtered pulses can be achieved by a combination 
of analogue delay lines, adders, and programmable pulse generators. 

[023] As shown in Figure 3, our invention can also apply a set of oscillators 
3 10 to the output of the filters 120 to shift the frequency, if necessary. The 
oscillators are of the form cos{Q\t + Frequency shifting enables an 
additional degree of freedom in the design, other than the delaying and 
weighting as shown in Figure 2. 

[024] An important part of our method is the determination of the pulse weights 
and delays. We determine the weight and delays on either fixed or time-varying 
criteria for the spectmm shape of the resulting signal 102. These criteria can 
stem from the FCC spectral mask (fixed), from the necessity to avoid 
interference to other users, which can be pre-defined or time-varying, or 
following an instantaneous or averaged determination of the emissions of users 
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in the current environment, or other criteria. In any case, these criteria are 
mapped onto an "instantaneous' spectral mask that has to be satisfied by the 
resulting signal 102. 

[025] Selecting Pulses Off-Line 

[026] The selection of the basis pulses 101 can be based on two criteria: 
implementation complexity, and approximation capability. 

[027] We prefer the basis pulses to be Gaussian in form. Gaussian pulses are 
relatively easy to generate, and their differentiations and integrations can be 
implemented with differentiation and integral circuits well known in the art. 

[028] The approximation can be performed with multi-resolution time- 
frequency analysis of signals using, for example, wavelets, and multi-resolution 
time-frequency approximation. 

[029] The signal of interest for the invention conforms to a desired spectral 
mask. Whereas a signal may be analyzed with respect an infinite number of 
possible bases in signal space, the choice of bases matters critically when one is 
restricted to a given number of bases available to approximate the spectrum of 
the desired signal. 

[030] A high-level description of the problem is given as follows. Assume there 
is a set of candidate of basis functions C. For example, the basis functions 
correspond to differentiations of a Gaussian pulse of various orders. In addition. 
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the basis functions can be orthogonalized, e.g., by a sequential Gram-Schmidt 
procedure, to facilitate the subsequent optimization. 

[031] A set of typical masks available for training purpose is denoted by 
S = {M(Q)} , and one particular set of chosen bases as denoted by The 
particular set is obtained by stacking the selected basis functions into a column 
vector. 

[032] We define a function / (^(0, S)to evaluate the fulfillment of the two 
considerations. For example, the function is a weighted combination of two 
parts: the first part models a cost of generating the pulse, and the second part 
models an approximation error or some other efficiency metric. This function 
can be expressed as the formulation: 

f{p{t), S) =af,{p{t)) + AiPiO, M (Q)), 

M (fi)€S 

where a and p are predetermined constants. This formulation is a combinatorial 
optimization problem, described below. 

[033] In comparison, the problem of optimal pulse locations is mainly for on- 
line applications, where the choices of bases are fewer, the implementation cost 
is fixed, and there is only one target mask, instead of a set of training pulses. 

[034] Formulation of the Problem 

[035] The individual basic pulses are denoted by p,(t), I = I,... L , and their 
Fourier transforms by i^OQ), / = 1,...Z . The set of shaping filter 120 is FIR, with 
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the impulse response being the sum of <5-functions placed at different delays, r, 
and weighted or scaled, 5, differently. Corresponding to Figure 2, we have 

1 = 0 

L Ml 

l=\ /=0 

L Ml 

°° 1=1 1=0 

We use the following notations: 

P(t) = [Pl(t-T,^) ... Px(t-T,^) ... P^it-T,^) ... Piit-TLMjf^ 

£p(OpitYdt, 

(s{t)Aty>= rs(tfdt=/Rs. 

J-OO 

[036] The elements of p(t) constitute a pool of bases pulses in the signal space. 

Then, the single user spectrum shaping problem can now be formulated as 
follows: 

max {s(t), s(t)), subject to \ S(jQ) f<M(Q), Vn e [-00,-1-00] , 

s 

where M(Q) is an upper-bound on the squared magnitude response regulated by 
FCC. 

[037] This is equivalent to a min-max formulation: 

I 5 ( >Q ) I' ^. r „ 

mm max — — — , subject to s R s = \ . 

s_ n6[-«. .«> ] M (Q ) — — 
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[038] In addition, structural constraints can be put on the parameter vector s 
due to implementation concems. For example, the number of non-zero elements 
in certain sub-vectors of the parameter vector s can be constrained. 

[039] In certain scenarios, the inner product of the desired 
signal (5 (O5 s (t)) can be replaced by the maximum of the signal energy 
within a certain frequency range, e.g., the "in-band" proportions of energy are 
maximized. This only changes the definition of the matrix R. 

[040] In certain scenarios, the constraints | S{jQ) f < M(Q), VQ e [-co,+co] 
are extended to include integral spectrum constraints, e.g., for "out-of-band" 
signals, which produces interference to other devices. These constraints are not 
as strict as the fixed envelop definition of the FCC mask. 

[041] The modifications corresponding to the system shown in Figure 3 follow 
after noting: 

fe-^^ cos(Q,r + (p,)Mt - T,)dt = i/?0(^^- Q,)K^^"-"' V^' + li?0(^^+ Q,)K'^"'"' V^' 
•^«> 2 2 

[042] Initializing Pulse Parameters 

[043] The min-max formulation as described above, or a robust 00-norm 
minimization is known to be a difficult problem. Several existing game- 
theoretic techniques rely on the existence of saddle points, which unfortunately 
are not satisfied in this case. 
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[044] Therefore, we minimize with an approximate 2-norm formulation 
instead. In other words, we replace the maximum over all frequency by the 
minimum of a weighted integral formulation, i.e., 

rl Si iQ) 1^ T 
w(Q)- — — — —dQ, subject to s Rs = l . 

We also note that: 
^ Af(Q) 



T 



[045] Provided that there is no structural constraint on s, the quadratic 
approximation leads to: 

min W s^, subject to s_ Rs^ - \ ^ 

for which the solution is an eigenvector corresponding to the smallest 
eigenvalue X of the matrix W-XR. When there are structural constraints on s 
and the non-zero positions have been fixed, e.g., as one elementary step in a 
combinatorial optimization procedure, the problem bears a similar form except 
that ^and R are now replaced by their corresponding principal sub-matrices. 
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[046] Optimization over Pulse Positions 

[047] One notable difference of our invention, compared with traditional least- 
squares FIR filter design, is that the pulse positions, i.e., delays of the non-zero 
elements in s are also optimized. 

[048] In prior art UWB systems, the random pulse-positive modulation relies 
essentially on a sparse arrangement of the pulses to support asynchronous 
transmissions, to provide better multi-path resolution, and to avoid catastrophic 
collisions. 

[049] The optimization over the pulse positions or delays is a combinatorial 
optimization problem, which determines an optimal subset selected from a 
given set of basis functions 

[050] As stated above, this problem is in similar form as the optimal choice of 
the initial basis functions, i.e., Pi{t - r^^) , from a larger set of candidate pulses. 
The second cost function f2{p{t),M{Q)) can be selected as an extreme value of 
the min-max, or the approximate min-integral formulation, 

[05 1] Here, we describe three selection procedures, with the first two being 
greedy. 

[052] Iterative Greedy Addition 

[053] We initialize with an empty set of basis functions. 
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[054] We search for a basis function from the candidate pool, which maximizes 
the cost function when added to the existing set of basis functions. 

[055] We select this basis function, and we repeat the search until the desired 
number of basis functions has been selected. 

[056] Iterative Greedy Removal 

[057] We initialize with the entire set of basis functions. 

[058] We search for a basis function from the candidate pool, which maximizes 
the cost function when removed from the existing set of basis functions. 

[059] We select this basis function, and we repeat the search until the desired 
number of basis ftinctions has been selected. 

[060] Branch-and-Bound 

[061] This approach assumes that p,{t - r„) have already been orthogonalized 

and normalized. Additionally, this approach uses a quadratic approximation 
rather than the min-max formulation. In this case, R = I, and the minimum 
generalized eigenvalue of matrix pair {W, R) reduces to the minimum 
eigenvalue of which is the cost function to be minimized. 
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[062] Let K denote the maximum number of non-zero elements. Let N denote 
the number of candidate pulses. Denote one arrangement of the pulse positions 
(delays) by a vector with ^ Ts, e.g., (1, 0, 0, 1, 0). 

[063] In the following, we constmct an enumeration tree of all the allowed 
combinations, where each node n is associated with a positioning vector with 
each element from the set {0, 1, *} where * means the position may be a 1 or 0. 
The enumeration tree is defined as follows: 

[064] The root of the tree is (*,*,. ..,*), referring to all possible selections. 

[065] The first layer of the tree enumerates all the vectors with one 0 and N - 
l*'s by considering the placement of a 0. 

[066] The children of any node with less than N - K zeros enumerates all the 

possibilities of putting another 0. To avoid duplications, ( , 0, *,...,*) is by 

convention associated with one possible father node ( ,*,*,...,*), where 

the first dotted part is the same as the current child vector. 

[067] At each visited node, the procedure maintains a lower bound and an 
upper bound of the cost function among all the descendant combinations. The 
procedure visits the enumeration tree with certain order, e.g., depth-first, or 
breadth-first. When visiting a certain node, the procedure updates the upper- 
bound and lower-bound of the cost function. If the lower-bound is greater than 
the minimum upper-bound among the visited nodes, the sub-tree, with its root 
being the current node, is pruned from further consideration. The bounds are 
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determined by Cauchy 's interlacing theorem of eigenvalues for symmetric 
matrices: 

Cauchy 5 Interlacing Theorem: 

Let Ar denote arXr submatrix of an X n real symmetric matrix A. The 
eigenvalues \ <X^<-'<X^ of A, and the eigenvalues /d^< ju^<'-< fi^of Ar satisfy 
the following relations: X. < < / = 1,. . .,r. 



The procedure can terminate after visiting all the nodes. In this case, the 
procedure reaches a global optimal selection of pulse delays. Altematively, the 
procedure can terminate when the cost function evaluated at a visited node is 
within £ of the current upper-bound of all the possible combinations. 



[068] Quadratic Optimization for Z = 1 



[069] For the special case of Z = 1 and r,. = / Ar, / = 1,...M, , a distinct quadratic 
approximation is pursued. In this formulation, the optimization over the pulse 
positions is directly solvable. 



[070] We introduce 

\p{m 

Then, 

„ax m= max 
n4-«^+«,] M{Q) G(Q) 
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Note, max ^ -1= max — — 

ne(-«,«o] G(Q) ne[-oo,-H»] G(Q) 



[071] Now, we select a substantially large constraint, i.e., s Rs = b^ such 
that the above equation is equal to the oo-norm. Although it does not make any 
difference in the original oo-norm formulation, it does affect the 2-norm 
approximation. With the 2-norm approximation, the problem becomes a FIR 
filter design problem with a least-squares formulation: 

min l^ll + s, ,e''^' + . . . + s.j^e''^'^' \ -G(fi) |^ t/Q, subject to /Rs = b with t oo 
[072] This equation can be approximated by 

A/, 



min V (5, - gi^„)^, subject to s^Rs = b with Z> t 00 



,=0 



where g„ is a real sequence that produces an amplitude response 
At 



[073] In particular, we match the non-zero positions to the taps of g„ with 
largest magnitudes. After, the positions of the taps have been fixed, the 

T 

solutions to the tap weights are computed as ^ " , w , g^^ where u„ is the 

principal component of R, and g are the vectors constructed from elements of g„ 
having largest magnitudes among those compl)dng with the structural 
constraints of s. 

[074] Note, the above two quadratic approximations are only two out of the 
many possible quadratic approximations. In essence, without FIR constraints 
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and the degree of freedom constraints, the optimal solution is an IIR filter, 
which uniformly matches the FCC spectral mask. 

[075] The various quadratic approximations can be viewed as minimizing some 
distance measure between the pursued FIR solution and the IIR solution. We 
prefer the above quadratic approximation mainly because the joint 
optimizations over the weights and delays are easily solvable due to its special 
stmcture. 

[076] Non-linear Optimization with Neural Networks 

[077] Figure 4 shows a feed-forward network with a differential soft- 
maximum, 410, and FIR frequency-responses at irregularly spaced frequencies 
420 to evaluate the function. In Figure 4, the weighting or scaling coefficients 
are s 401, the pulse positions or delays are r 402, 0/ represents | (Q, , s, r) |^ , e"^ 

431 represents the exponential operator, and !P^41 1 represents the output of the 
soft maximization, i.e., the sum of the exponentials 43 1. 

[079] We initialize our solution with the above quadratic approximation, and 
further exploit non-linear optimization techniques to gradually refine the 
solutions. The description in this section refers to the original min-max 
formulation described above. 

[080] A back-propagated (BP) multi-layer perceptron (MLP) possesses 
adaptive leaming abilities to estimate sampled functions, represent these 
samples, encode structural knowledge, and inference inputs to outputs via 
association. Its main strength lies in its substantially large number of hidden 
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units, and thus, a large number of interconnections. The MLP neural networks 
enhance the ability to learn and generalize from training data. 

[081] We describe the MLP optimization for the special case L = I. 
Optimization of the scaling coefficients, with fixed positions in the general case 
L> I, follows similarly. Because Z = 1, all quantities that depends on L are 
simplified to shortened notations with the dependence removed. 

[082] In addition, we re-parameterize the problem so that s/, 5jt represent 
the scaling or weighting coefficients at the pulse locations (delays) r,,. . .r^ . 

Therefore, we define L = [^i>"''^kT . The optimization over pulse positions 
for Z, = 1 can be treated by MLP, for Z > 1 . We refer again to the three 
combinatorial optimization approaches discussed earlier, 

[083] In order to put the current problem into the general framework of MLP, 
we uniformly quantize the frequency range to arrive at 

|0o41 max ■ — — — « max - — ^^^^ — '-LJ^ 
n6[-«.,+«>] M(Q) .•6(0,1.... Af-i) M (Q i) 

[085] rthrh 



Then, we replace the max-ftinction with a differentiable soft-max. Given a 
function fix) >0,xeS ,fora sufficiently large positive number a, we have 
the following soft-max approximation: 

[086] max f(x) « —In Y 



17 



MH.5145 
Molisch et al. 



Then, we replace the max-flinction with a differentiable soft-max. Given a 
function f{x) > g 5 , for a sufficiently large positive number a, we have 
the following soft-max approximation: 

max /(jc) « —In Y e^'^^'K 

[087] With these two simplifications, the problem becomes 

T 

min y/ z^), subject to s_ R {r_)s_ = 1 , 

with R{t) being the principal submatrix of R corresponding to delays r and 

/=: 0 

[088] Figure 5 shows the feed-backward network 500 for adaptive optimization 
with back-propagation, where 501 and 502 are the derivatives defined by the 
equations in Fig. 5. 

[089] While it is theoretically possible to simultaneously adjust both the 
positions or delays r 40 land the weighting or scaling coefficients s 402, 
practically we prefer to decouple their tuning by adopting a conditional 
maximization approach, i.e., optimizing one with the other fixed. In addition, 
this decoupling may be justified by the different nature of the two parameter 
sets. In preferred embodiment, with the allowed positions quantized, the slight 

change I^^T+M^r is always normalized to "hop" to the nearest valid 
quantization point on the multi-dimensional grid. 
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[090] Typically, numerical non-linear optimizations can only be assured to 
arrive at a local optima rather than a global one. Simulated annealing can be 
used to avoid local optima. One approach to escape from possible local minima 
is to begin the optimization with several randomly distributed initial solutions, 
and select the best solution among the different trial paths. 

[091] We note that this is only one possible neural network, which again is just 
one of the methods for implementing non-linear optimization. 

[092] Although the invention has been described by way of examples of 
preferred embodiments, it is to be understood that various other adaptations and 
modifications may be made within the spirit and scope of the invention. 
Therefore, it is the object of the appended claims to cover all such variations 
and modifications as come within the true spirit and scope of the invention. 
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